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The influence of (non-magnetic and magnetic) impurities on the transition temper- 
ature of a d-wave superconductor is studied anew within the framework of BCS theory. 
Pairing interaction decreases linearly with the impurity concentration. Accordingly 
suppression is proportional to the (potential or exchange) scattering rate, 1/r, due to 
impurities. The initial slope versus 1/r is found to depend on the superconductor con- 
trary to Abrikosov-Gor'kov type theory. Near the critical impurity concentration Tc drops 
abruptly to zero. Because the potential scattering rate is generally much larger than the 
exchange scattering rate, magnetic impurities will also act as non-magnetic impurities as 
far as the Tc decrease is concerned. The implication for the impurity doping effect in high 
Tc superconductors is also discussed. 
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1. Introduction 



Recently there has been much attention on the d-wave pairing state in connection 
with high-temperature superconductors. ^'^'^ The impurity effect on the d-wave pairing is 
particularly interesting because it may give a clue to the symmetry of the superconducting 
state of the high cuprates. There are already many theoretical works on this issue. ^"""^^ 
For examples, the impurity effects on the penetration depth, ^ the superfluid density,^ the 
transition temperature,^ the quasiparticle states,^ the density of states, ^'■'^^"■'^^ the nuclear 
spin relaxation rate,^ the spin susceptibility,^ and the infrared conductivity,^ have been 
investigated. 

The above theoretical works were essentially based on the Abrikosov-Gorkov's (AG) 
pair breaking theory^'^ for the s-wave superconductors. However, Kim and Overhauser^^ 
proposed a new theory with different predictions: (i) The initial slope of decrease 
depends on the superconductor and is not the universal constant suggested by AG. (ii) 
The Tc reduction by exchange scattering is partially suppressed by potential scattering 
when the mean free path is smaller than the coherence length, which has been confirmed 
in several experiments. Kim-*^® also showed that the failure of AG theory comes from 
the intrinsic pairing problem in Gor'kov's formalism, which can be cured by incorporating 
the pairing constraint. The purpose of this paper is to reconsider the Tc suppression due to 
(non-magnetic and magnetic) impurities in a d-wave superconductor within the framework 
of BCS theory. 

It is shown that the pairing interaction decreases linearly with the (non-magnetic 
or magnetic) impurity concentration. The ratio of the average correlation length, = 
0.18^, to the mean free path, ^0/^5 determines the weakening of the pairing interaction. 
Consequently Tc suppression and the initial slope depend strongly on the superconductor 
contrary to Abrikosov-Gor'kov type theory. For the magnetic impurities, because the 
cross section of the potential scattering is generally much larger than that of the exchange 

2 



scattering, the decrease is also determined by the potential scattering of the magnetic 
impurities. ^^'^^ 

The interpretation of the high Tc data^*^'^^ by the present study is not clear, though. 
The observed Tc decrease due to impurity scattering is much slower than the theory pre- 
dicts. The discrepancy seems to come from the neglect of the strong Coulomb interaction 
effect in calculating the impurity response of the superconducting state. In other words, 
the impurity responses are different in Fermi liquids and strongly correlated systems. Note 
that it was pointed out^^ the Tc decrease caused by impurity doping and ion-beam-induced 
damage in high Tg superconductors is related with the proximity to a metal-insulator tran- 
sition. 

2. D-Wave Superconductor 

For a d-wave superconductor, the pairing interaction ^, for the plane states is taken 
to be^^ 

= J e'^'-''^-^>(r)rfV = -5V2 ■ k'f - 1], (1) 

where % is the unit vector parallel to k. Substituting Eq. (1) into the BCS gap equation, 
one finds 

= E • k'f - ^]^^tanh^, (2) 

fc' ^' 

where 

= (3) 
and is the electron energy. Among the possible solutions, we consider 

A,- = Ao(fc2-fcJ). (4) 

This solution has the same symmetry property as d^^-y^ = Ao{coskx — cosky) which 
is believed to describe the gap structure of the cuprate high Tc superconductors. Upon 
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substitution in Eq. (2), one obtains the Tg equation 

Te = 1.13eee-l/^°^^ (5) 

where ec is the cutoff energy and No is the electronic density of states at the Fermi level. 

In the presence of impurities, the scattered states i/jn may be expanded in terms of 
plane waves, such as^^ 

V;„ = ^e^^-'"<^|n>. (6) 
k 

Now the pairing interaction Vnn' between scattered basis pairs (V'niV'n) ^ind (V'n')V'n') is 
given by 

Vnn'= [ I dTidT2rn'iri)rn'ir2)V{\Vi-r2\)Mr2)M^l)- (7) 



Here ipn denotes the time-reversed state of ijjn- From Eqs. (1), (6) and (7) we can calculate 

Vnn'- 

3. Non-Magnetic Impurity Eflfect 

The non- magnetic impurities will be examined first. The scattering potential from 
the impurities is given 

U{r) = Y^u5{r-Ri). (8) 

i 

{Rj} is the impurity sites. We consider the impurity effect using the 2-nd order perturba- 

— * 

tion theory. Then, the scattered state which carries the label ka is 



i,q ^ k-\-q 



where A^^ is the normalizing factor. The time-reversed degenerate partner of Eq. (9) is 



= + E e'^-^'e-'^''+^-'']p. (10) 

i,q k k+q 



The matrix elements which cause Cooper pairing are between scattered basis pairs 
ii^ka^^-kp) (V'fc/Q,' ^-fc'/j)-"^^'^"^ Each basis pair is a 2 x 2 Slater determinant. The 
matrix element Vj^ p between the two scattered basis determinants is 

Vk,k' =< 'Dk'i^i,r2)\V{\n-f2\)\V^{rur2) >, (11) 

where 



i^-kf^i^i) ^-kp{^2) 



(12) 



Upon employing Eqs. (1), (11), and (12), we find 



Vk,k' = -5^2 \m ■ k'f - mM,. (13) 



The normalizing factor N'i is given 



Nr = \ nr, (14) 



where \ Wj^\ is the relative probability contained in the virtual spherical waves surrounding 
the impurities (compared to the plane wave part). As in s-wave superconductor case, in 
computing we cut off the radial integral at i? = 3.5^o/2, because the pair-correlation 

amplitude falls exponentially as exp{—r /3.5^o) near Tg.^"^ The average correlation length 
is defined by 

Co = 0.18^. (15) 
With Eqs. (13), (14) and (15), one readily finds 



1 

I 

where £ is the mean free path. 

Therefore the pairing interaction is reduced: 



Vk,k' = -^^^ ^[3(^^0^-l][l + ^]-^. (17) 



Notice that in dilute limit the reduction is proportional to the ratio of the average correla- 
tion length to the mean free path, ^oji and the pairing interaction decreases linearly with 
the impurity concentration. The Tc equation is now, 

Te = 1.13eee-V^°^^[l+^l-^ (18) 

Figure 1 shows versus l/r for T^o = 40K and 80K respectively. Tco denotes the 
transition temperature without impurities. We used ec = BOOK. For a metal with vp = 
2 X lO^cm/sec, the superconductivity is completely suppressed when the mean free paths 
are about lOOOA and 350A for Tco = 40K and 80K, respectively. As in magnetic impurity 
effect on s-wave superconductors,-^^ we find a sudden drop of Tc near the critical impurity 
concentration. This effect was much weakened by potential scattering in s-wave case. In 
this case, this effect may be real. Note that Porto and Parpia^^ found a sudden drop of Tc 
in p-wave superfluid He^ caused by aerogel. On the other hand, the inelastic scattering 
may decrease the effect in high Tc superconductors. 

The change in Tc relative to Tco can be calculated to first order in the impurity 
concentration. From Eqs. (15) - (18) 

A^c«^-^, (19) 

where r is the collision time due to impurities. The factor 1/ A shows that the initial slope 
(versus l/r) depends on the superconductor and, consequently, is not a universal constant. 
Weak superconductors lose their Tc more rapidly than strong ones. 

4. Magnetic Impurity Effect 

The magnetic impurities give rise to both potential and magnetic scattering. For the 
potential scattering, the result of the previous section is applicable. In this section, we 
consider the effect of magnetic scattering. The magnetic interaction between a conduction 
electron at r and a magnetic solute (having spin S), located at Rj, is 

Hm = Js-SiVoSir-Ki), (20) 
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where s = |cr and Vq is the atomic volume. The three components of a are the Pauh 
matrices. As in Sec. 3, we consider the magnetic impurity effect up to the second order of 
J. Then, the scattered state which carries the label, ka, is 

= ^fj;^-he''■'a + J2 e'^'+'^-{Wj^^^ + w'^.a)], (21) 



where. 



and. 



= '^^^^ E -nx,e^^--^^ (22) 



'^L-=i^?^E-«^--'"'- (23) 



k k + q j 

Xj and are the polar and azimuthal angles of the spin Sj at R-,-, and S = -^f S{S + 1). 
The degenerate partner of Eq. (21) is: 

= -A/fe^"^ [^"''"^/^ + E e-'^'^^-^\WlM - Wy.(5)]. (24) 

q 

In this case, we pair and 

Accordingly, the matrix element ^, is given 

^M' =< ^fe'(^"'i>^"'2)|nk"'i -^"'2|)|2^fe(ri,f2) >, (25) 

where 



*_fe;3(ri) *_fc^(r2) 



(26) 



Using Eqs. (1) and (26), we find 

Vk,k' = -5^2 Im ■ k'f - 1]AA|AA|. (27) 

Note that Eq. (27) is the same form as Eq. (13) in Sec. II. As a result, the Tc reduction 
caused by the magnetic impurities can be calculated by the same formulas, Eqs. (17) - 
(19). The only difference is using ig and Tg instead of £ and r. £g is the mean free path for 
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the spin disorder scattering and Tg is the spin-disorder scattering time. However, because 
the cross section for the potential scattering of the magnetic impurities is generaUy much 
larger than that for exchange scattering, we only need to consider the potential scattering 
due to magnetic impurities as far as Tc change is concerned. 

5. Discussions 

This study has been done in the framework of the BCS theory. But the same re- 
sult can be obtained from the Gor'kov's Green's function theory if we impose a pairing 
constraint on the self-consistency equation. Near the d-wave self-consistency equation 
with degenerate pairing constraint is given by^^ 

A*{T,l) = Vir-\)Tj2 [ [ dr'dl'{GM{r',T,-cu)GN{l',l,oj)}''A*ir',l'). (28) 

w 

The superscript P denotes the pairing constraint and Gn is the normal state Green's 
function in the presence of ordinary impurities. Using the relation between A„ and A* (r, 1), 



A„=/ drd\A%r,\)i;nmn{r), (29) 




it can be shown that Eq. (28) is the nothing but another form of the BCS gap equation. We 
can also show that the physical constraint of the Anomalous Green's function F^(r, 1, a;), 
(i.e.), 

-imp — T-p -imp 



Ft(r,I,c.)""- = Ft(r-l,a;) , (30) 

gives rise to the degenerate pairing constraint. The superscript ~~^'^p means an average 
over the impurity positions. 

In high Tc superconductors, the impurity doping and ion-beam induced damage^^ 
suppress strongly Tc. But the Tc reduction is not fast enough to be explained by this 
study. The experimental data show that Tc reduction is closely related with the proximity 
to a metal-insulator transition caused by the impurity doping and the ion-beam-induced 
damage.^°'^^'^^ We believe that the impurity response of the strongly correlated systems 



like high cuprates may be understood only when we consider the impurity scattering 
and the strong correlation on equal footing. Because both non-magnetic and magnetic 
impurities suppress Tc almost equally, anisotropic pairing may be plausible. 

6. Conclusion 

In conclusion, we considered the effects of the (non-magnetic and magnetic) impuri- 
ties on Tc in a d-wave superconductor. The pairing interaction decreases linearly with the 
impurity concentration. The initial slope of the Tc suppression depends on the supercon- 
ductor and therefore is not a universal constant. Because the cross section of the potential 
scattering is much larger than that of the exchange scattering, Tc suppression is deter- 
mined entirely by the potential scattering of the magnetic impurities. We also discussed 
the implications of this study for the impurity doping effect in high Tc superconductors. 
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Figure Captions 

Fig. 1 Variation of Tc with impurity concentration (measured in terms of the scattering 
rate, ^) for T^o = 40i^ and 80K, respectively. The cutoff energy ec is 500 K. 
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